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Abstract 

In this paper, we use the Weisskopf-Wigner theory to study the entangle- 
ment in the state of the free-space radiation field produced from vacuum 
due to atomic decay. We show how bipartite entanglement is shared be- 
tween different partitions of the radiation modes. We investigate the role 
played by the size of the partitions and their detuning with the decaying 
atom. The dynamics of the atom-field entanglement during the atomic 
decay is also briefiy discussed. Prom this dynamics, we assert that such 
entanglement is the physical quantity that fix the statistical atomic decay 
time. 
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Entanglement plays a central role in quantum information science where it 
is seen as an important physical resource for information processing beyond 
the achievable with classical correlations [T||2||3] . Investigations in quantum 
phase transitions [^115] and statistical mechanics have been realized from the 
point of view of entanglement [6] . In this approach, concepts like randomness, 
ensemble-averaging or time-averaging are not required. Instead, thermaliza- 
tion results from entanglement between system and environment. Connections 
between matter and quantum information theory have also been discussed 
[7118] . Such studies have pointed out that entanglement seems to be important 
in other areas of physics besides pure quantum information. It is exactly in 
this scope that this work is introduced. 

Considering entanglement as a legitimate physical quantity, this paper is in- 
tended to study a fundamental process namely the atomic spontaneous emis- 
sion. It lies at the core of matter-radiation interaction. The successful descrip- 
tion of atomic decay in free space is one of the remarkable achievements of 
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the quantum theory of radiation pilOllllj . Why does an excited atom decay? 
It is clear that an isolated atom would never decay from one excited state to 
another with lower energy because both are eigenstates of the system Hamil- 
tonian. There must be some physical system to couple to the atom in order 
drive the electronic transition. This external agent is the free space electro- 
magnetic field whose zero-point energy fluctuations are able to cause the atom 
to decay. In the language of quantum information theory, the entanglement 
between atom and field is then the responsible for the atomic decay. In what 
follows, entanglement is studied in the spontaneous emission phenomenon. 

We are particularly interested in the entanglement properties of quantum 
fields. In this paper, the quantized field is the bosonic free-space continu- 
ous electromagnetic field. We will analyze the entanglement properties of this 
field after atomic decay. Entanglement in discretized bosonic fields have al- 
ready been studied. In particular, the entanglement properties of the ground 
state and thermal states of this system is studied in detail in [12]. Discrete 
versions of real free Klein-Gordon fields have also been studied from the point 
of view of entanglement ^LS}. In this study, the relation between entanglement, 
entropy, and area for a specific harmonic lattice whose continuum limit lead 
to the Hamiltonian of the real Klein-Gordon field is analyzed in great detail. 
These papers triggered many others studies that also investigated entangle- 
ment in discretized quantum fields P^15II16|[T7] . 

The starting point of the present work is the Weisskopf-Wigner theory of spon- 
taneous emission which is now briefly presented pilU|llj . In the rotating wave 
approximation, a two-level atom interacts with the free-space electromagnetic 
field according to the interaction picture Hamiltonian [10] 



H = hY^[gl{ro)a^a^e^^--"^^' + H.c], (1) 

k 

where u is the angular frequency of the atomic transition (excited state \a) 
and ground state \b)), is the position of the atom, at is the annihilation 
operator for the field mode {k} (angular frequency Uk), and 



^7k(ro)=^?ke-^'^-'-°, (2) 
where 



with V a quantization volume, eo the electric permittivity of free space, pab 
the dipole moment for the atomic transition, and ek the polarization vector 
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of the mode {k}. It will be assumed that initially the atom is in the excited 
state \a) and the field modes are in the vacuum |0) = |0,0, ...). According to 
([T]) the system evolved state will be 

= Ca{t)\a, 0) + E Cfe,k(t)|&, Ik, {0}), (4) 

k 

where |lk, {0}) represents the field state with one photon in the mode {k} 
and the rest in the vacuum, and 

Ca(t)=e-^*/^ (5) 
c,k(t)=^k(ro) (,^_,) + ,r/2' 
with r being the free-space atomic decay constant which is giving by 

1 AuJ^fp^„^. 



In order to obtain the above equations, it was considered that the intensity 
of the light associated with the emitted radiation is very centered about the 
atomic frequency u. This is the essence of the Weisskopf-Wigner theory. In 
this theory, the free space modes act as an immediate response reservoir, i.e., 
the atomic spontaneous emission is seen as a Markovian process. 

Now, the entanglement content in the field state after spontaneous decay of 
the atom is studied in detail. This state is denoted |7o) and it is obtained from 
01]) by assuming t ^ 



l7o) = E , ,°^p;J lk,{0}). (8) 

It is worth noticing that the state I70) is, from the point of view of quantum 
information science, a member of an important class of multipartite entangled 
states called generalized W states [18]. However, we must take care when using 
this state. Although state ([8]) is presented as a discrete summation over k, any 
kind of calculation using it is to be done transforming it to an integral, i.e. an 
continuum of modes. 

The state ([8]) represents all modes of the free space radiation field, and it is a 
superposition of the different possibilities of distributing one photon (emitted 
by the atom) between the infinity of modes. Consequently, this is an entangled 
multipartite state whose bipartite entanglement between different partitions 
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of radiation modes is now going to be investigated. There are many ways of 
partitioning the free space modes in two partitions. We think it is physically 
appealing to choose one partition formed by a central mode with frequency z/q 
and modes distributed in the interval (z/q — e, z^q + e) (let us call it partition A), 
and the other partition formed by the rest (partition B). This is an interesting 
physical choice since it allows us to check the effect of having z/q either near 
or far from resonance with the decaying two-level atom (frequency separation 
u), and to check the importance of the size of the partitions via the parameter 
e. 

Since |7o) is a pure state, the appropriate entanglement measure between 
partitions A and B of the system is the entropy of entanglement E = S{pa), 
where S{pa) = — tr[p^ log2(pA)] is the von-Neumann entropy with the reduced 
state Pa = tT^ b[pab]- It must be emphasized that the entropy of entanglement 
is a entanglement monotone only if the global state is pure. Even though the 
pure field state |7o) is achieved only in the limit t ^ F^^, we will see later on 
this paper that our results are still approximately valid for finite times. This 
broadens the applicability of our work. The reduced state for the partition A 
can be obtained from ([8]) by tracing out modes in partition B. One finds 



PA=Ebri{o}A)u{o}i+ E p-p:iik^>{o})(ik.,{o}i, (9) 

kj km 1 kn 

where kj refers to a wave vector of some mode in the partition B, km (kn) to 
some mode in partition A, \{0}a) to vacuum states of modes in partition A, 
|lk^(^), {0}) means one photon in mode km(n) of partition A and vacuum for 
the rest of the modes in that partition, and 



(z/,,-o;) + zr/2- ^''^ 

The only non-zero eigenvalues of pa are Ai = X^kj bjP ^2 = Skn IPn]"^, 
where kj refers to partition B and kn to partition A. As mentioned before, 
the final results must be obtained by passing to the continuum. In spherical 
coordinates we have [10] 
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For consistency with the Weisskopf-Wigner used in the derivation of the state 
([8]), we should again consider that varies httle around z/^ = uj, what allows 
us now to replace u'^ by uj^ in the above integrals as well as to extend the lower 
integration limit of the first integral to — oo llOj. Making this approximations 
one obtains 



1 arctan 

TT 



Now, we sum the modes referring to partition A 



— arctan 

TT 



(13) 



Pab 



(14) 



Again, we replace by (but leave the integration limits unaltered) to 
obtain 



E 1^*" !^ ~^ ~ arctan 

It is important to look into normalization of ([H]) because we must end up with 
a physical state after performing the approximations. In fact, the normaliza- 
tion has been conserved since the sum of f|T3l) with f|T5|) is equal to one for 
any values of e, co and Uq. With f|T3|) and f|T5l) . one can now easily obtain the 
entropy of entanglement S = — Z]?=i -^j study the bipartite entan- 

glement between partitions A and B. In general, two special features of the 
entanglement in the field modes should be highlighted, namely its dependence 
upon the size of the partitions and upon the detuning between the central 
frequency of partition A and the atom A = Ug — u. From now on, we will use 
the dimensionless quantities e = e/F and A = A/F. 

In Fig.([T]), we show how the entanglement varies with the size of partition 
A. We can see that the effect of increasing the size is initially to increase 
the entanglement between both partitions. It physically means that more and 
more entangled modes are shared by the partitions. One would expect the 
entanglement to saturate because as e increases it comes to a point where 
the number of modes in each partition optimizes the bipartite entanglement. 
If e continues to increase, the situation returns to be unbalanced with less 
and less entangled modes shared by the partitions. In the limit e — > cx), there 



UJ] 



— arctan 

TT 



(15) 



are no modes left in partition B, and the entanglement between the partitions 
naturally goes to zero. In Fig.([T]), one can also see that if the central frequency 
i/q moves from resonance with the atom (increase of |A|), it takes more and 
more modes (larger values of e) to achieve maximum entanglement as expected. 
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Fig. 1. Entanglement between partitions A and B as a function of the size of parti- 
tion A. The curves corresponds to different values of A: 0.0 (solid), 2.0 (dash-dot), 
4.0 (dot), 8.0 (dash). 

We have just seen that the entanglement between two partitions of the free- 
space modes after atomic decay considerably depends on the detuning between 
the atom and the central mode of partition A. Such dependence is show in 
Fig. ([2]). First, one can notice that if the size of the frequency band e is small, 
the maximum entanglement is not achieved for any detuning. This is in full 
accordance with all that has been already discussed here about the role played 
by e. On the other hand, if e is sufficiently high, there is always a value of A 
which allows the system to achieve the maximum entanglement. It is physi- 
cally expected that if one is to consider a partition whose central frequency is 
far from resonance with the atom, the modes will almost not be affected by 
the atomic decay. In fact, one can see in Fig. ([2]) that in the limit A ±oo the 
entanglement goes to zero. We should make it clear that A can not actually 
be considered too big since we are in the scope of nonrelativistic quantum 
mechanics in which high energy interactions are not properly accounted. In- 
deed, even the Hamiltonian ([T]) obtained in the dipole approximation would 
not be valid in such regime. In spite of that, the Weisskopf-Wigner theory 
used here is very accurate in the optical-microwave domain and fully explain 
the main features of atomic spontaneous emission in free space. An attempt 
to go beyond the Weisskopf-Wigner theory is presented in [19]. It is also im- 
portant to remark that momentum entanglement of the atom and the field is 
not included in our treatment. Such problem has been solved in [201I2T] where 
the authors conclude that such entanglement is very small indeed. 

The Weisskopf-Wigner theory then allows one to learn about the entanglement 
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Fig. 2. Entanglement between partitions A and i? as a function of the detuning 
between the atom and the central mode of partition A. The curves corresponds to 
different values of e: 0.2 (solid), 0.5 (dash-dot), 5.0 (dot), 9.0 (dash). 

in the field modes after atomic spontaneous emission. Since this theory also 
gives the time evolution of the global state comprising the atom and the field, 
one can go a step further and obtain the dynamics of entanglement between 
these two subsystems as well. From Eqs.(jl]), ([5]), and ([7]), we may easily obtain 
the reduced density operator for the atom (tracing out the modes), and it is 
given by 



Pat(t)=e I a) (a 



k 



Just like before, we now pass to a continuum of modes, and the result after 
integration is 



p.t{t) = e-^'\a){a\ + {l-e~^')\h){hl (17) 

which coincides with the atomic density operator obtained in the master equa- 
tion formalism in the Born-Markov approximations [22]. In situations where 
the Markov approximation breaks, the Weisskopf-Wigner approach might be 
quite useful for the the study open quantum systems. An interesting problem 
where the Weisskopf-Wigner theory applies is inhibition of atomic spontaneous 
emission in photonic cristals which is known to be a non-Markovian process 
[23] . Another special feature of the Weisskopf-Wigner approach is that it al- 
lows us to study the state of the reservoir as well. For atomic spontaneous 
emission in free space, the study of the reservoir was performed in the first 
part of this paper. 

The reduced system has now just two states, and the evaluation of the entropy 
of entanglement is straightforward. Its time evolution is shown in Fig. ([3]). 
The maximum of entanglement takes place at tm = ln(2)/r. Consequently, 
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Fig. 3. Time evolution of the entanglement between the atom and the electromag- 
netic modes of free space. 

we can reinterpret the maximum of entanglement as the physical quantity 
that fix the time needed for half population in an atomic ensemble to decay 
spontaneously (this time is called half life in nuclear physics). At t^, and 
for the particular problem treated here, the atom is in a Bell state with the 
reservoir [S'(rtjn) = 1]. 

At first sight, the results presented here seem to come into conflict with what 
is expected in the Born-Markovian approximations which is known to be valid 
for atomic spontaneous emission in free space. For example, the Born approx- 
imation states that the atom-reservoir state becomes only slightly different 
from a perfect separable state [22|. Consequently, the entanglement shown in 
Fig-© should not be so high. In fact, if we were to evaluate the entanglement 
between the atom and each mode or with a small group of modes, we would 
see that there is really just a small deviation from product states. However, 
when we consider the whole bunch of free space modes, this bipartite entan- 
glement may become surprisingly high. Consequently, even though it is right 
to say that the atom-reservoir state is almost separable due to the fact the 
reservoir is a very large system (not much affected by the atom), it is exactly 
because it is that big that the small entanglement between each mode and 
the atom leads to high entanglement when all modes are considered, as shown 
in Fig. ([3]). We think this a central result in our paper. The validity of the 
factorization assumption is quantitatively studied in a perturbative approach 
in [21]. 

When we think of Markov approximation in this problem, i.e. the indepen- 
dence of the future behavior of Pat(^) on its past history, it is usual to justify 
it on the following reasonable physical grounds [22] • If the reservoir is a large 
quantum system, it is expected that it is mantained in thermal equilibrium, 
so that it is not supposed to preserve the changes induced by the atom at 
previous times. Since we are dealing with the zero temperature case, the ther- 
mal equilibrium state reads the vacuum state. Now, it is our findings in the 
first part of this paper, i.e. the study of the entanglement properties of the 
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field state after spontaneous emission, tliat seem to come into conflict with the 
Markov approximation. Again, the size of the reservoir is key for understand- 
ing it. The field state after atomic decay ([H]) contains just one photon shared 
by an infinity of radiation modes. We physically expect that this state is in fact 
very close to the vacuum state at the beginning of the atomic decay process 
(except for the fact that it is an entangled state). If this is true, the Markov 
approximation is still valid because of the very small deviation from the ther- 
mal equilibrium. A naive calculation of the fidelity F = Tr(|7o) (7o|pvacuum)5 
would lead to F = 0. This is not relevant for the kind of problem we are 
dealing with (infinitely many subsystems). For example, we expect the states 
|0, 1) to be "more different" from |0, 0) than |0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) is from 
|0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0). Basically, our expectation is based on the number 
of subsystems in the vacuum state. The more subsystems in the vacuum, the 
more such a state is close to the global vacuum state. However, the fidelity is 
zero for both cases in this example. Again, it is physically more appealing to 
concentrate our attention in one partition instead of the global state |7o). In 
Fig.(jl]), we show how the state (Q of the partition A compares to the vacuum 
state by evaluating F(e, A) = TrA(pAPvacuum) in the continuum limit. We can 
see that there is a vast region in which the fidelity is equal to one with respect 
to the vacuum. Although it is not a formal proof that the field state stays 
at equilibrium, one can easily see that depending on the cut or partition, the 
resulting si ' ' ' " ' 




Fig. 4. Fidelity between the state of partition A and the vacuum as an function of 
the size of the partition and the detuning between its central mode and the atomic 
frequency. 

We would like to return to the discussion of the validity of our results for 
finite times. The entanglement properties of the field state studied in this 
paper were obtained under the assumption that the overall state of the modes 
was pure. We now argue that this high purity is still approximately true even 
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for finite times. We first recall that the distillable entanglement En [25] is 
lower bounded by entropies as follows [3] 

Ed{pab) < max[^(pB) - S{pab), 0], (18) 

where, in our case, A and B are partitions of the field modes. In the spe- 
cial instance of a overall pure field state [S{pab) = 0], the above inequality 
becomes an equality and entanglement is quantified by the entropy of the re- 
duced state, as we did in this work. In order for our results to be valid (at 
least qualitatively), we must make sure that the entropy of the field modes 
S{pab) = is in fact approximately zero. From the Araki-Lieb inequality [25] 
S{pa.t,AB) < \S{pab) — S{ps.t)\^ and the fact that the atom-field state (jl]) re- 
mains pure all the time, we see that S{pab) = 5'(pat)- Fig- ([3]) clearly shows 
that entropy of the atom (consequently the entropy of the modes AB) goes to 
zero quite rapidly. For example, S{pab) = 10~^ is achieved for Tt ^ 7. As one 
can see, our results remain essentially valid even for finite times, and this fact 
may facilitates the experimental investigation of the entanglement properties 
discussed here. 

To summarize, we have studied bipartite entanglement in the atomic sponta- 
neous emission phenomenon. For the field state after atomic decay, we have 
found that very detuned modes weakly entangle with the rest, and that the 
size of the frequency band of each partition strongly changes the entangle- 
ment. We have carefully analyzed the behavior of entanglement when such 
parameters change, and also physically justified all features. It is important 
to remark that the present work represents a study about fundamentals of 
matter-field interactions from the point of view of quantum information, and 
it does not necessarily aim at applications. Actually, applications of the en- 
tanglement in the free space radiation field state after simple atomic decay as 
considered here seems unlikely. Of course, atomic decay in free space could be 
used to entangle distant atoms as shown in P7f25] . However, what we have 
just studied here is entanglement in the free space modes. This makes coherent 
control a challenging task when compared to atomic entanglement. 

On the other hand, we believe that a scheme to experimentally access the field 
entanglement induced by atomic decay could make use of photonic crystals. 
Spontaneous emission in this environment has been studied in [23], and since 
this system is much more controllable than the electromagnetic fields in free 
space, we feel that this could be an alternative to study photonic entanglement 
due to atomic decay. It is hoped that the present work may be a starting point 
for addressing this question (theoretical treatment of entanglement properties 
of atomic decay in photonic crystals) as well as other interesting related ques- 
tions such as the application of the techniques used in [19j to study atomic 
decay. 
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